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Abstract -We analyze the condensation phase transitions in out-of-equilibrium complex networks 
in a unifying framework which includes the nonlinear model and the fitness model as its appro- 
priate limits. We show a novel phase structure which depends on both the fitness parameter and 
the nonlinear exponent. The occurrence of the condensation phase transitions in the dynamical 
evolution of the network is demonstrated by using Bianconi-Barabasi method. We find that the 
nonlinear and the fitness preferential attachment mechanisms play important roles in formation 
of an interesting phase structure. 
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' Introduction. — Condensation phenomena emerge 
in various physical contexts, to name a few, the well- 
known Bosc-Einstein condensation (BEC) in dilute atomic 
gases [IH3], jamming in traffic flow [HIS], wealth conden- 
sation in macroeconomics [B], and condensation in zero- 
r,ange process (ZRP, see e.g., recent review [7] and ref- 
erences therein). Since the pioneered research on com- 
plex networks [5HTB]. in the last decade, condensation phe- 
nomena, i.e., condensation of links (or edges) in complex 
networks has also been widely discussed [TTHHS] . In the 
Context of complex networks, the condensation phase cor- 
r,esponds to the situation that a single node captures a 
ihacroscopic finite fraction of total links/edges. It has 
been found that condensation phenomena can occur in 
both growing and non-growing complex networks. The 
condensation phase transitions occurring in non-growing 
networks [20l[52l[23l[25ti28| are formally equivalent to that 
in balls- in-boxes model [35] , and hence has been well stud- 
ied to a large extent via methods of equilibrium statistical 
mechanics. While for growing complex networks, the ap- 
pearance of the condensation phase transitions during the 
dynamical evolution of the network are particularly inter- 
esting due to its out-of-equilibrium characteristics. 

The tasks in this paper are two folds: first, we merge 
two important models on this regard, the growing net- 
work with nonlinear preferential attachment (we will refer 
to "nonlinear model" [TH| from now on), and the fitness 



model [37] into a unifying framework-the nonlinear fitness 
model. We then argue that the condensation phase tran- 
sition appearing in fitness model and that in nonlinear 
model stems from different mechanisms; Second, particu- 
larly interesting, we reveal a novel phase structure in the 
model and this may increase our understanding on the 
non-equilibrium phase transitions in dynamical evolution 
of complex networks. 

The nonlinear model is defined as follows. At each time 
step t, the newly-added node created m directed links to 
ones of the earlier existing nodes with fc-link, according to 
a probability, say, 11, that is proportional to some "connec- 
tion kernel" k'^ ,Ii (xk^ . Here the exponent 7 > refiects 
the tendency of preferential linking to a popular node and 
hence controls the preferential attachment. In Rcf. [T5] . 
P. L. Krapivsky et al. had discussed the cases of different 
choices on exponent 7, i.e., 7 = 1, 7 < 1 and 7 > 1, for 
growing complex networks with connection probability. 
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They proved that the number of nodes with k links, N^, 
follows a power law distribution in the case that 7 closes to 
unity. While in the case of 7 < 1, the distribution shows a 
stretched exponential form. For 7 > 1, that is so called the 
super-linear case, the model exhibits a condensation phase 
transition. Especially when 7 > 2, there exists a limiting 
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situation where the most connected node links to almost 
all the other nodes in the network, this corresponds to a 
"winner-takes-all" phenomenon. In this case, the degree 
of the most connected node follows kmax t. 

A similar condensation phase, or so-called "winner- 
takes-all" phase, also appears in fitness model of the com- 
plex networks [TUII^IST] . In fitness model of growing net- 
works, a fitness parameter, ry^, which represents an inter- 
nal superiority of the i-th node, is introduced and chosen 
randomly from some distribution. As a result, the well- 
known Barabasi-Albert scale-free network [5] is general- 
ized to the Bianconi-Barabasi (B-B) fitness model of com- 
plex network. One may assign the i-th node the fitness 
rji according to its "energy level" Si which satisfy some 
distribution g{e) through the relation 
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where /? can be identified as inverse temperature, i.e., /3 = 
1/T. In this growing complex network with fitness, the 
connection probability II-'' that a new node connects one 
of its m links to an existing node i at each time step t is 
defined by 



(3) 



where ki is the degree (the number of links occupied by one 
node) of node i. One may introduce a partition function 
Zt as 



(4) 



This model can be solved in a mean-field approximation 
and the condensation phase transition process is described 
by the "chemical potential" /i formally defined as follows. 



lim — 

t-i-oo mt 
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where Zt is the partition function Zt averaged over some 
normalized distribution 5(e) |19j. 

The chemical potential is determined by the self- 
consistent equation 



= / deg{e)n{e) = 1, 



(6) 



where n[e) is the occupation number, i.e., the number 
of links attached by the preferential attachment mecha- 
nism to nodes with "energy" e. Very interestingly, it was 
proved [T^] that n{e) = l/(e^(^~^^ — 1), this is nothing 
but Bose-Einstein (BE) statistics. When /i < 0, the net- 
work is in so-called "fit-get-rich" (PGR) phase. While in 
the thermodynamic limit t ^ 00 if I{/3,0) < 1, i.e., the 
self-consistent equation ([6]) has no solution, and hence a 
BEG phase transition occurs. There exists a critical tem- 
perature Tc = l/Pc such that /(/3,0) < 1 for T < Tc- 
This condensation phase transition was demonstrated as 
well by numerical simulations in Ref. [T^ . 



Both the fitness model and nonlinear model of the com- 
plex network experience the condensation phase transi- 
tion during their dynamical evolutions, two questions are 
naturally arisen: whether or not the condensation phase 
transitions in these models have some common underlying 
relationship and the role that different preferen- 

tial attachment mechanism in cither model plays during 
the phase transition of the network. In order to answer 
such questions, and as well, to explore the possible phase 
structure in out-of-equilibrium evolution in complex net- 
works, we propose the nonlinear fitness model of complex 
networks which includes the nonlinear model and the fit- 
ness model as its appropriate limits, as we will show in the 
following. 

The nonlinear fitness growing network. It is 

clear that the preferential attachment mechanism controls 
the network topology in the nonlinear network and the 
fitness parameter rji does in the network with fitness. To 
bridge the divergence between the two models and to give 
a general description of the phase structure, we employ a 
new connection probability, ^l'^ , as follows. 
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This new probability includes the original connection 
probability 11^ in the nonlinear model and II-'^ in the fitness 
model, respectively, as its proper limits: in the limit 7 — >^ 1 
one goes back to the preferential attachment probability 
of the fitness model, and in the limit 77 — >■ 1, one recovers 
the nonlinear model of growing network. This attachment 
mechanism now addresses the dynamical evolution of the 
network. In order to identify the condensation phase tran- 
sition of this nonlinear fitness model of network, similar 
to the Bianconi-Barabasi method, (see Eqn. ([5])), we for- 
mally define the chemical potential fji [T5] . 
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where, /3 is again the inverse temperature, m is the 
number of newly added links at each time step t, and 
Zt = J2j=i^~^'^^^"'i^j^^T^j) is partition function in 
our nonlinear fitness model, t — > cxd plays the role of ther- 
modynamic limit. 

The corresponding rate equation is. 



dt ~ ^ ' 



(9) 



Note that in 7 — 1 limit, the rate equation of B-B fit- 
ness model is restored and one naturally expects a BEG 
phase transition at low enough temperature [in]. However, 
in general, the chemical potential /Li(7, T) is a function of 
both temperature T (or the fitness 77) and the nonlinear 
exponent 7, i.e., the phase structure of this nonlinear fit- 
ness network is controlled by these two parameters, in- 
stead only one parameter in either model (T in the orig- 
inal B-B fitness model or 7 in the nonlinear one). These 
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parameters both affect the phase transition during the dy- 
namic evolution of the network. This fact leads to a more 
complicated behavior of fj. in current model than that in 
the fitness or the nonlinear model alone. Here we apply 
the method of rate equation of degree, used in the B-B 
fitness model |19) . rather than that of connectivity distri- 
bution, used in the nonlinear model fl8| . This is due to the 
fact that the fitness parameter rji is different for each node 
i, i.e., 77 is a local, not a global parameter. This suggests 
that the fitness parameter 77 and the exponent parameter 
7 which respectively controls the preferential attachment 
mechanism are different in nature - 77^ reflects an "inner" 
property of the ith node, while 7 controls the global evo- 
lution of the network. This important difference between 
the B-B fitness model and the nonlinear model leads to 
the different dynamical evolution result of the complex 
networks. 

Numerical simulations and discussions. Simi- 
lar to the original B-B fitness model, we identify the non- 
condensation-condensation phase transition by the change 
of sign of the chemical potential /i(7, T), i.e., when /j, ex- 
periences a change from negative value to positive value 
within some regime, that change implies the critical point 
of the corresponding phase transition. In addition, in prin- 
ciple, the phase transition occurs in the thermodynamic 
limits t — > 00. However, in general situations, taking such 
thermodynamical limit is not realistic, one has to resort to 
the numerical simulations. We numerically compute the 
chemical potential /it according to Equ.®. 

Throughout our simulations we fix the number of links 
per node (per time step) m = 2. For simplicity, we take en- 
ergy level distribution g(e) = Ce with normalization con- 
stant C = 2, and the total number of time steps t = 10^, 
average over 100 runs. The main numerical simulation re- 
sults are plotted in a three dimensional (3d) figure (see 
Fig. [1]), in which three axes of the figure are exponent 7, 
temperature T and chemical potential /i, respectively. In 
this "phase diagram", one can see that the role that the 
temperature T or exponent 7 plays in the formation of 
such a phase structure during the dynamical evolution of 
the network. For instance, for 7 = 1, when one lowers 
the temperature, a transition from FGR phase (high tem- 
perature phase) to BEC phase (low temperature phase) 
occurs [in] as what we expect. Similarly, if one fixes the 
temperature, say T = 2.0, a phase transition would also 
occur when the exponent 7 is varied (between 7 = and 
7 = 1.5 in Fig.[T|). However, one should keep in mind that 
the current phase transition is different with that in non- 
linear model because no fitness (temperature) effect was 
taken into consideration in the latter model. On the other 
hand, in both cases, the signal of the phase transition lies 
on the change of sign of the chemical potential. Note that 
we do not take absolute value of the chemical potential, 
which is different with Ref. [H]. We see that in Fig. [T] the 
chemical potential /i experiences a change from positive 
values to negative ones when the exponent 7 and/or the 
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Fig. 1: (Color online) Three dimensional illustration of the 
phase structure in the nonlinear fitness complex network (Top), 
and its projection on the 7 — T plane (Bottom) . In top panel, 
it shows the chemical potential /i as a function of temperature 
T and the nonlinear exponent 7, as defined in Eqn. ((§]). In this 
simulation we take the number of links per node (per time step) 
m = 2, and energy level distribution g(e) — 2e (normalized), 
and the total number of time steps t = 10'^, average over 100 
runs. 



temperature T vary in some regime. Now the whole phase 
structure is richer than that in the B-B fitness model and 
in the nonlinear model, respectively, e.g., the BEC phase 
and the FGR phase in the B-B fitness model now are only 
parts of the new phase diagram in 7 = 1 limit. 

The situation is more clear when we plot the chemical 
potential /i as a function of exponent 7 in Fig. [5J and 
of temperature T in Fig. [31 respectively. These figures 
display the different "cross section" views of 3d "phase 
diagram" of Fig. [TJ Some new features of the phase tran- 
sition can be read out from these plots. In Fig.[2l we show 
the chemical potential ^ versus the exponent 7 with tem- 
perature T fixed. It can be seen from the figure that as 
temperature increases, the critical value of exponent for 
phase transition 7c decreases, e.g., for T = 0.2, jc ^ 1-2, 
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Fig. 2: (Color online) The chemical potential /i(7, T) as a func- 
tion of the exponent 7 but with temperature fixed, in a linear- 
linear scale. The four curves shown in the panel correspond 
to T = 0.2 (open triangles, solid line), T = 0.8 (open circles, 
dashed line), T — 2.0 (open squares, dotted line), T = 5.0 
(open diamonds, dot-dashed line), respectively. All symbols 
are connected by lines for eye guidance. 

while for T = 5.0, jc ^ 0.6, roughly. This tells us that 
the transition of phase occurs with weaker tendency of 
preferential linking at relatively higher temperature. 

In Fig. |3] we plot the chemical potential fi versus the 
temperature T with exponent 7 fixed in a linear-log scale. 
For data points with 7 = 0.7 (open circles, dashed line) 
and 7 = 1 (open squares, dotted line), the condensation 
phase transition is very similar to what happens in Fig. [21 
the critical temperature Tc of the phase transition de- 
creases as the exponent 7 increases, e.g., for 7 = 0.7, 
Tc 3.0, while for 7 = 1, Tc 0.8. The latter case 
recovers the results of B-B BEG phase transition, as it 
should. Fig.|3]also displays some new, interesting aspects: 
the values of the chemical potential with 7 = 0.3 (open 
triangles, solid line) are all positive and those with 7 = 2 
are all negative. This implies that the phase transition 
disappears eventually as 7 reaches some specific critical 
value, and hence results in a single phase structure - ei- 
ther a FGR phase (with negative chemical potential), or 
a condensation phase (with positive chemical potential). 
However, one can not always use this as an identification 
of transition between the condensation phase and the non- 
condensation phase, this point is illustrated in Fig. 

We plot the occupation ratio, kmax/'mt, of the most con- 
nected node, in Fig. At the first sight from Fig. SI it 
seems that the network is in a condensation phase at low 
enough temperature T < Tc {Tq is the critical tempera- 
ture), and in a non-condensation phase at high tempera- 
ture, T > Tc- But note that as the exponent 7 increases, 
the condensation phase appears in even higher tempera- 
ture region. Especially when 7 = 2, the network enters 
a condensation phase regardless of the value of temper- 



Fig. 3: (Color online) Similar to Fig. [2l the chemical poten- 
tial ^(7, T) as a function of temperature but with exponent 7 
fixed, in a linear-log scale. The four curves shown in the panel 
correspond to 7 = 0.3 (open triangles, solid line), 7 — 0.7 
(open circles, dashed line), 7 = 1.0 (open squares, dotted line), 
7 — 2.0 (open diamonds, dot-dashed line), respectively. Note 
that for 7 = 0.3, the values of /i-s are all positive, while for 
7 = 2.0, the values of /i-s are all negative. 



ature. This supports the conclusion of nonlinear grow- 
ing network in Ref. |18| . However, in addition to the ap- 
pearance of the condensation phase in original B-B fitness 
model and the nonlinear model with 7 > 1, we conclude 
from Fig. H] that they are still different phases because the 
chemical potential is positive > 0) for the B-B fitness 
model at r < Tc, while negative {fi < 0) for the non- 
linear model. The competitive relationship between e~^^^ 
and k'^ in partition function play important roles in the 
dynamical evolution of corresponding complex networks 
and hence the formation of such a novel phase structure. 
When the former factor e~^^^ is dominant, the network 
experiences a condensation phase transition if T < Tc, 
and its chemical potential is positive, > 0; while the 
prevailing role of the factor k'' leads to the opposite side: 
the chemical potential of the condensation phase is always 
negative, n < 0. It is also clear that, from above figure, 
for /i > and 7=1, the network with fitness preferen- 
tial attachment is condensed on the node with the lowest 
energy level (or the fittest node) only if T < Tc in the 
thermodynamic limit; for /i < 0, the network experiences 
a "winner takes all" phenomenon when 7 is large enough 
(e.g., 7 > 1.7). Due to aforementioned competitive rela- 
tionship between two factors e~^^'- and k'', the chemical 
potential of the non-condensation phase is changed, for 
instance, see the curve with 7 = 0.3 in Fig. |3l though 
one has an overall positive fi at high temperature, no con- 
densation phase appears. The different phase structures 
identified by the chemical potential /i stem from different 
dominant factors of preferential attachment mechanism, 
they are distinct in nature. In addition, as is well known, 
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Fig. 4: (Color online) The occupation ratios of the most con- 
nected node, kmax/mt, plotted as a function of temperature T, 
and of different nonlinear exponent 7. The five curves in the 
panel, from bottom to top, represent 7 = 0.3 (open triangles, 
solid line), 7 = 0.7 (open circles, dashed line), 7 = 1 (open 
squares, dotted line), 7 = 1.3 (stars, solid line), 7 = 2 (open 
diamonds, dot-dashed line), respectively. Note that the data 
points with 7 = 1 are consistent with B-B's results (see the 
corresponding curve of Fig. [S] in Ref . [19] ) . 



a significantly different manifestation of BEC and FGR 
pfiases is tliat the fraction of the total number of links 
connected to the most connected node tends to zero in 
thermodynamic limit for the latter and tends to a finite 
value for the former. In current combined model. Fig. 2] 
shows that there exists a similar situation for ^ > when 
the factor -q dominates the attachment mechanism, while 
when the factor 7 dominates the attachment, the ratio 
tends to zero in thermodynamic limit for /i < 0, which 
displays a FGR phase. However, note that if 7 is large 
enough, this occupation ratio can also approach to a finite 
value and a condensation phase appears even for ^ < 0. 

A major difference between the condensations of fitness 
model and non-linear model is that the condensation in 
the former occurs on the node with the lowest energy level 
(or the highest fitness), instead in the latter model it only 
occurs on the first node of the network. However, topo- 
logically it is hardly able to tell whether or not a graph 
has reached condensation due to the non-linear preferen- 
tial attachment or the BEC due to the fitness, since in 
both condensation phases, the important topological pa- 
rameter, i.e., the cluster coefficient (CC) approaches to 
the same limit. If one trace the node with the lowest 
energy level, then this node is more likely to be the most 
connected node in low temperature region, since with large 
parameter /? = l/T, the fitness parameter 77 dominates the 
dynamical evolution in combined attachment mechanism. 
As temperature increases, the probability of the node with 
the lowest energy level being the most connected node is 
getting smaller due to the impact of fitness 77, as a result. 



Fig. 5: (Color online) The occupation ratios of the most con- 
nected node, kmax/nit, and that of the node with the lowest 
energy level, k^/mt, plotted as a function of temperature T, 
and of different nonlinear exponent 7's. Data with hollow sym- 
bols and solid lines correspond to kmax/nnt, represent 7 = 0.7 
(open circles), 7 = 1 (open squares), 7 = 2 (open diamonds), 
respectively. Data with solid symbols and non-solid lines cor- 
respond to k^/mt, represent 7 = 0.7 (solid circles, dashed line), 
7=1 (solid squares, dot-dashed line), 7 = 2 (solid diamonds, 
dotted line), respectively. 



in high temperature region, the factor k'^ dominates the 
attachment. This is shown in Fig. [5j in which the oc- 
cupation ratios of the node with the lowest energy level, 
k^/mt, is plotted; as well as that of the most connected 
node, kmax/nit, as function of temperature. The former 
exhibits evidently larger fiuctuations compared to the oc- 
cupation ratio of the most connected node kmax/nT-t, since 
the node with the lowest energy level is not necessary the 
most connected node. As a result, in high temperature 
region, in which fc'*' dominates the attachment (77 — >■ 1), 
fluctuations are relatively small, since as a global param- 
eter, 7 does not depend on the local fitness of a node, 
such that the earliest presented node could be the most 
connected node for large 7, regardless the local fitness for 
T > Tc (Tc is the critical temperature); in low tempera- 
ture region, (T < Tc), fluctuations of the occupation ratio 
of a node with the lowest energy level suggests that the 
emergence of condensation comes from the competition of 
two sides - the local fittest node (with fitness attachment) 
and the earliest presented node (with super-linear prefer- 
ential attachment), and this competition between makes 
fluctuations larger and larger (e.g., see 7 = 2 in Fig. [S]). 
This demonstrates how both attachment factors t] and 7 
control the evolution of the network simultaneously. 

Conclusions. — In conclusion, in this paper we an- 
alyze the condensation phase transitions in a unifying 
framework which includes both the nonlinear model and 
the fitness model as its appropriate limit. The goodness of 
this new framework is to, on one hand, bridge the differ- 
ences between the nonlinear model and the fitness model; 
On the other hand, allow us to identify the different roles 
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played by different factors in current model. Under ap- 
propriate limits (e.g., 7 — > 1 or 77 — ^ 1), the typical phase 
structures of the original B-B fitness model and the non- 
linear model are recovered respectively. 

Through the numerical simulations of this nonlinear fit- 
ness model, we show a novel 3d phase diagram. We employ 
the Bianconi-Barabasi method and numerically compute 
the chemical potential to identify the critical points {Tc 
and 7c) of the phase transitions of the networks. Our re- 
sults and analysis directly answer the questions that asked 
in the beginning, i.e., the condensation phase transitions 
in both the B-B fitness model and the nonlinear model are 
in fact distinct in nature. As well, the preferential attach- 
ment mechanism in our nonlinear fitness model depends 
on two factors: k'' and e~^'' (or rji). Both factors affect 
the phase structure of the network and hence the network 
topologies during the evolution of the network. We reveal 
that the competitive relationship between these two fac- 
tors leads to the condensation phase transitions, and hence 
transitions between different network topologies. We hope 
the current study may increase our understanding to the 
non-equilibrium critical phenomena, particularly the con- 
densation phase transitions in the evolution of complex 
networks. 
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